Constraints on the dark energy equation of state from recent supernova data 
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Measurements suggest that our universe has a substantial dark energy component. The most 
recent data on type la supernovae give a dark energy density which is in good agreement with other 
measurements if the dark energy is assumed to be a cosmological constant. Here we examine to 
what extent that data can put constraints on a more general equation of state for the dark energy. 

PACS numbers: 98.80 -k, 98.80.Es 



1. INTRODUCTION 



It is generally believed, based on the analysis of type la supernovae red shift data that the expansion of 

the universe is accelerating and that this acceleration is either caused by a "dark" energy or is the effect of extra 
dimensions. In the case of dark energy this energy could be a simple cosmological constant term in Einstein's 
equations or it could be a time dependent energy density. In any case the pressure of the dark energy must be 
negative to account for the acceleration. What we would really like to know, to determine the nature of past 
and future acceleration, is the equation of state of the dark component, X, as a bmction of the scale parameter 
R{t). This equation of state is defined as 

w{z) ^ ^, (1) 
Px 

where px is the pressure, px is the energy density, and z is given by the ratio of R{t) to the scale today, Rq, as 

1 + z = Ro/R{t). (2) 

For a cosmological constant w{z) ~ —1 for all z. If the acceleration is caused by extra dimensions it can also 
be parameterized as an effective w{z). 

There now exists a new data set of measurements of the luminosity distances for type la supernovae, which 
range in distance from z = 0.01 to z = 1.755 This data includes some new events as well as a reanalysis of 
the old events. There were about 194 old events 0, Q some of which have been dropped and most of the rest 
assigned new errors. The new set consists of 157 "gold" events or 186 gold plus "silver" events. 

The gold data set gives joint confidence levels in the matter -dark energy plane with w = —1 that are more 
restrictive than the earlier data and in good agreement with other constraints This is shown in Fig.^ Thus, 
perhaps, it is not too early to ask whether this data can distinguish between simple models for w(z). To investigate 
this question we use a few models such as w{z) = wq + wi z, w{z) = wq and w{z) = wq + wi z/(l + z) to 
calculate a luminosity distance. We then minimize the of the difference between the measured and calculated 
luminosity distances in an attempt to determine the constant parameters wq and wi. WMAP and other 
experiments have accurately determined the matter density and curvature, defined as 



rig rjQ 



to be close to = 0.3 and flk = 0. po is the density of matter today, Hq is the Hubble constant today, and 
k ~ — 1,0, +1. Nevertheless, because we would like to know to what extent the supernova data agrees with 
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FIG. 1: (Color online) In the left panel, the 68.3% (red), 90% (dark blue), 95.4% (blue) and 99.7% (green) confidence 
contours for the gold data fit with a cosmological constant are shown. The cross indicates the best fit, and the solid circle 
the best fit to a flat cosmology. The right panel shows the 68.3%, 90% and the 95.4% confidence contours for the data in 
Ref. Q. 

these other experiments, we will sometimes include f^M and/or f2fe in our set of parameters to be fit despite the 
fact that this exacerbates the degeneracy problem that arises because the luminosity distance depends only on 
the sum of the matter, curvature and dark energy contributions. 

In an earlier paper we explored the sensitivity of the data in Refs. to simple assumptions about 
the dark energy density. Recently, many other papers have appeared which use the data of Refs. 

to analyze numerous models of the dark energy. A paper which is close in spirit to this work and which we had 
the advantage of reading before we started this work is Ref. While we disagree with some of the results 
and conclusions, there is still much to like about this reference; many more models were studied and physical 
motivation is given for some of the models. We, on the other hand, have chosen to consider only the simplest 
cases to remain as model-independent as possible. 

Our results are not encouraging. The errors in the fitted parameters are often more than 100%, in some cases 
much more. We estimate the error in itself is more than 10%. Still some consistencies can be seen and some 
patterns discerned. And while the data is not manifestly in agreement with other experiments (e.g., = 0.3 
and f2fe = 0) it is not in disagreement. One (negative) result that seems very clear is that the absolute value of 
X^ cannot be used, at this time (i. c., from this data set), to discriminate between models for w[z). 



2. CALCULATION OF THE LUMINOSITY DISTANCE AND ^ 



The luminosity distance is given by 

dL{.z) = 

where H{z) is the Hubble parameter 

HI 



c(l + z) 



sm 
1 

I I sinh 



(5) 



(6) 
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nx = 1 - - ilk, (7) 
and f{z) is given by the equation of state 



w{x)) . (8) 



In Eq. (|SJ) sin[ ] is used for flk < (fc = +1), sinh[ ] is used for ilk > (fc = —1). and the unmodified square 
bracket is used for fife = 0. 

Many references give models where f{z) is not obviously related to an equation of state. However, from Eq. JHJ 
we see that /(O) = 1, f{z) > 0, and, given an explicit expression for f{z), 

1 f'(z) 

w{z) - -1 +3(1 (9) 

where the prime denotes a derivative with respect to z. This equation can also be used to express the effects of 
extra dimensions as an effective 'w{z). 

Below we use the dimensionless luminosity distance 

D,{z) = (10) 
c 

and the quantity that is compared to the data is 5 logj^p Dl{z) + M where M is a constant offset between the 
data and the theoretical expression. The comparison is made by calculating where 

^.2 - ^ logio-Pr''(^0 - - 5 \og-yQDL{zi)f ^^^^ 
The sum is over the N = 157 events contained in the gold set 0] and cr,; are the experimental errors in 

5iogioi^r(^o- 

Dl{z), through the equation of state w(z), Eq. J^l, depends on some set of parameters, ai, 02, ... , a„, and we 
minimize by varying these parameters, s-lso depends on M but to minimize with respect to M is trivial, 

N 

(12) 



where N and D are given by 



^ ^ J2 5 1ogiogr(^.) ^5 log,, Dl{z,) ^ ^^g^ 



1=1 

N 

D = 



^ 1 

1=1 * 

Sometimes we use Eq. (|12() for A/, sometimes we include M as one of the parameters Ui to be fit, and wc have 
checked that we get the same result. 

The data indicates that the universe was decelerating in the past and it is of interest to know when that 
changed. The acceleration is given by 

(15) 



1 dluH^ 



1, (16) 



2 d ln(l + z) 

and, given from Eq. ij^Jl, and values for the parameters, it is straightforward to find where 9(2,) = 0. 
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3. RESULTS 

We consider three simple, rather ad hoc, expressions for the equation of state 

z 

w{z) ^ wq + wi— — , (17) 
i + 2 

w{z) = Wo, (18) 
w{z) — wq + wiz , (19) 
where Wq and Wi arc constants. These expressions allow us to evaluate the integral in Eq. © and avoid having 
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TABLE I: Fit of supernova la data to w = wo + wi — wiy, y = {1 + z)^^ (first 3 rows), w = wo (rows IV, V and VI), 
w = Wo + wiz (rows VII, VIII, and IX) and some special cases (last three rows). The errors are the 68.3% confidence 
level numbers. The boxes without errors mean that parameter was held fixed. If a parameter was varied no restrictions 
were placed on its possible values except for Qm > and ^Ix > 0. The last column uses the values of the parameters and 
their errors to find the z value where the universe changes from decelerating to accelerating. 

an integral within an integral. We minimize with respect to wq, wi (if relevant), flm, (or 0.x), or some 
subset of these; sometimes we hold Vlk fixed at zero or Qm fixed at 0.3 or wq fixed at —1 or some combination of 
these conditions. If a parameter is allowed to vary we impose no prior conditions on the values it may assume 
except for > and fix > 0. 

It is not straightforward to minimize because the minimum is very shallow and very broad. We use two 
independent numerical codes to cross check, one is a slig ht g eneralization of programs (MRQMIN) given in 
Numerical Recipes the other is a code called STEPIT |13| which is available on the web. 
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Our results are given in TableQ] The first 3 rows use Eq. IjlTfl . with cither nothing, ilk, or flM and fife fixed. 
The next 6 rows do the same for Eq. (|18|l and then for Eq. (|19() . The last 3 rows give some special cases. For 
each model we give the minimum value and the values of the parameters which give that value. The errors 
on the parameters arc the 68.3% confidence level uncertainties given for a particular parameter by finding the 
variation in that parameter which results in an increase of Xahti by 1 when the other parameters are allowed 
to vary freely. Plots of how varies to the 1 a error in each parameter are shown in Fig.El for Models I (red 
curves)and VII (blue curves). If the parameters in these models are restricted to the case of a flat cosmology, 
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FIG. 2: (Color online) The variation of the parameters Q,m, fix, wq and wi is shown for the cases ■w{z) 
wi z/{l + z) (red) and ■w{z) — wo + wi z (blue). The horizonal dashed lines correspond to Ax^ ~ 1. 



Wo + 



fife = 0, the resulting 1 a errors of Models II and VIII are shown in Fig.O 

To try to get an estimate of the uncertainty in xi/m itself, and to check the stability of the numbers in Tabled 
we have done approximately 120 bootstrap simulations. These simulations consist of constructing new data sets 
by randomly selecting 157 events from the original data set. We do enough such simulations that each event 
gets used about the same number of times. We then calculate the average xi/m ^^r the simulations and, to find 



6 



176 



175 



X 174 



173 



172 



, , , , 1 , , , , 1 , , , , 1 , , , < 1 , , , , 1 , , , , 




\. 011=0.475 / 




^^^^ X^=173.6 / 




n„=0.459 


.... 1 ... . 


X^=173.5 

.... 1 .... 1 .... 1 ... . 



0.0 



0.1 



0.2 



0.3 
fin 



176 



175 



0.4 0.5 



176 



0.6 



175 



174 



173 



173 



X 174 



173 



172 



, , , , 1 


, , , , 1 , , , , 1 , , , , 1 






\ W(,=-2.44 / / 






\\ x''=173.6 // 






Wo = -2.4B 




.... 1 


X^ = 173.5 

.... 1 .... 1 .... 1 





.... 1 .... 1 .... 1 .... 1 .... 1 




Wi = 1.74 / / ■ 




^^=5^X^=173.6 / / - 




Wi=3.75 




X^=173.5 

1 .... 1 .... 1 .... 1 ... . 



-20 



-15 



-10 



-5 



10 



FIG. 3: (Color online) Same as Fig|5|except that f2i 



an error for xiiin j take the standard deviation of the set of Xm in ■ We also find the average and the standard 
deviation for each of the parameters. These results are given in TablellTI The values of the parameters and the 
errors in the parameters are very consistent with those given in TablcD In fact for the cases where only one or 
two parameters are fit they are almost identical. The values of for the individual simulations vary from 110 
to ~ 230 with a standard deviation of about 20. Thus, from these simulations, wc conclude that the values in 
Tableware robust and that Xuin is uncertain by 10 — 15%. 



CONCLUSIONS 



If we take the cosmological constant case {wq = — 1, wi = 0) and fit or fi^ we get the preferred answer of 
0.3 or (see Models X and XI in Tabled . Or if wc hold Hm ~ 0.3 and ^lk — and fit wq we get the preferred 
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TABLE II: The averages of approximately 120 bootstrap simulations for the models of Table I. The errors are the standard 
deviations. 

answer of —1 (Model VI). However, if we try to fit even two of the three wq, ^m, as in Models V or XII, the 
fitted values are not the preferred ones although the errors are large enough to make the parameters consistent 
with the desired values. This can be seen in some detail in the wq — rij\/ confidence contours for Model V, shown 
in Fig. 21 Model XII, shown in more detail in Fig.Q was our original motivation. 

Our Models III and IX are models 9 and 7 of Ref. ^3 ■ That paper used the 2003 data set 0, ; when we 
ran with that data set, and included the extra error due to the peculiar velocity as in |lOj |. we were able to find 
a smaller value for xifin both of these cases; perhaps we used a different expression for the pecular velocity 
error. In any case a smaller XMin ^'^^ very important in view of the large uncertainty in x\iin indicated by 
the simulations. What is important, in view of this uncertainty, is that small relative values of XMm '^ot 
indicate which models are a better fit. Also the ability to find the lowest is essential to determining the errors 
on the parameters; if the lowest is not found the errors will appear to be smaller than they actually are. 
We found much larger errors for the parameters than those given in Ref. |lOl | and this seems not to be due to 
uncertainity in handling the pecular velocity error. For example our result for Model III is wq = —1.74 ± 0.39, 
1^1 = 5.04 ± 2.25. 

If we include structure in w{z), as in Models I, II, III, which use w{z) = wq + wi z/{l + z) or VII, VIII, IX, 
which use w{z) = + z, then the errors get very large and it is hard to draw any conclusions even if fij\/ 
and f2fc are fixed. What does seem clear, from the similarity of the results for Models I and VII, II and VIII, 
and III and IX is that low z events are still dominating, [wi of Models VII, VIII, and IX should be compared 
with wq + w\ of Models I, II, and III.) The factor of (1 + z) in the denominator of Eq. lfT7|l is effectively just 1. 
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FIG. 4: (Color online) The 68.3% (red), 90% (dark blue), 95.4% (blue) and 99.7% (green) confidence contours for tlie 
case Qx = 1 — Sljv/ and px = wopx (model V) are shown. The cross indicates the best fit, and the solid circle is the fit 
requiring SIa/ = 0.3. The dashed line separates accelerating (to the left) from decelerating (to the right) cosmologies. 

We have checked that this is also true for the gold plus silver data set. Thus comparing models for the equation 
of state of the dark energy will remain something of a mug's game until there exists substantially more data at 
higher values of z. 
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